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1. Introduction 

Recently thermal field theory in the light front (LF) frame introduced by Dirac [|J] has gained 
quite some attention. The most important application in this framework is the phase diagram of 
strongly coupled systems like, e.g., the quark gluon plasma. Todays perception of the QCD phase 
diagram is due to by lattice QCD computations. However, these calculations are limited to the 
region T < fx (T temperature, jU quark chemical potential) due to the complex action at large 
chemical potential. In turn, this results in the well known sign problem of the Monte-Carlo simu- 
lation method. The generic Monte-Carlo sign problem is at least as hard to compute as problems 
in the complexity class NP (class of non-deterministic polynomial problems) and every problem in 
NP is reducible to the sign problem in polynomial time (i.e. the generic Monte-Carlo sign problem 
is NP-hard) [Q] and therefore we argue that it is worth looking for alternative ways to determine the 
QCD phase diagram. 

In the following we investigate LF quantization to compute thermodynamical quantities. The 
first attempt to use results of light front quantization, i.e. the invariant mass spectrum and the 
wave functions of the theory, for applications in thermodynamics has been given in |@|. How- 
ever, the conclusions of Ref. [|3j] are rather confusing since a second order phase transition in 
one-dimensional QED has been conjectured. The limiting cases of non-interacting fermions on 
one hand and the free boson gas on the other have not been considered. Certain classes of su- 
persymmetric models and four-dimensional pure gluonic QCD [Q] have been investigated and 
thermodynamical properties computed. Analytical calculations in LF thermal field theory have 
been performed for different models. These perturbative computations have been done using a 
statistical operator familiar from the more traditional instant form approach. It was possible to 
reproduce known results like thermal masses in scalar field theory and properties of the Nambu- 
Jona-Lasino model [||, ^]. A notation of the general light cone (GLC) frame, which compromises 
between instant and front form coordinates, was introduced [||] and in the following it was pointed 
out that the canonical quantization in the GLC frame [^] is essentially analogous to ordinary light 
cone quantization. However, the advantages of light cone quantized thermal field theory stemming 
from technical simplifications in perturbative computations like the simple pole structure of the 
propagator have been hardly exploited, see e.g. JTO)] - 

A non-perturbative approach to light cone quantized field theories is given by discrete light 
cone quantization (DLCQ) [[Tl|]. Discrete light cone quantization is a finite box quantization of 
Hamiltonian field theory supplemented by boundary conditions for the fields and cuts the Fock 
space into finite-dimensional sectors of equal resolution K = j^P + , where L is the box length. Mass 
spectra and LF wave functions of low lying states which are independent of the box length have 
been numerically computed for one-dimensional or dimensionally reduced systems via DLCQ. 
Higher dimensional systems are usually treated by the transverse lattice approach which replaces 
two spatial dimensions by a lattice and the remaining two by DLCQ. The problem of renormaliza- 
tion in Hamiltonian field theory and therefore the construction of effective light cone Hamiltonians 
remains to be solved and hampers application of light cone quantization to non-perturbative quan- 
tum field theory in 3+1 dimensions. 

In this proceeding we carefully reconsider questions in QEDi + i as raised in [Qj. However 
arrive at mostly different conclusions concerning the conjectured phase transition. Some of our 
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results have been recently given in Ref. [12]. 



2. Light Front Thermodynamics 

Starting from considerations in Ref. [ |l3| ] the statistical operator on the light front can be written 
in the following form 

1 



apparently different from the 'naive' light cone version pic ~ ex P(j3^~) that resembles the non- 
relativistic form. The partition function is given by 2£ = Trp which is the central quantity when 
one wants to compute thermodynamical and statistical properties. When evaluating the partition 
function in DLCQ one introduces the harmonic resolution K and the light cone Hamiltonian H 
through 

, 2k L L M 2 

p + = T K ' r = (2 " 2) 

Here K is dimensionless, diagonal in the DLCQ basis and used as a measure of the discrete ap- 
proximation. The light cone Hamiltonian has dimension mass squared and is the dynamical part in 
( |Q| ) since it is a non-diagonal matrix of increasing size in K. Inserting (2.2) the partition function 
reads 

^ ( r, i )=Trexp{§(^,/ + A«*§)}, (2.3, 

where 'Tr' means summing over all resolutions K and all corresponding (decoupled) Fock space 
sectors, and / is the identity matrix. The mass matrix Mk of course is different for different K- 



sectors. Note that the volume appears explicitly in (2.3) in contrast to the suggestion for i2° in [g]. 
This is due to the consistent approach based on eq. (2. 1). Mo is the mass of the lightest state in the 
continuum limit, that means we normalize the smallest eigenvalue of Mk to one for K — > oo. In a 
numerical computation we fix the volume (in units of the continuum estimate of the lowest mass) 
at the beginning and extrapolate for K to infinity. This calculation has to be performed for several 
values of L to safely determine the expected linear dependence 

Q. = -Tln2f = aL + p, (2.4) 

where Q. is the thermodynamical potential. We emphasize that one has to pick a strict order of 



limits in L and K, first take K — > oo followed by L — » oo. Computing (2.3) in practice means expo- 
nentiating large matrices and summing the diagonal elements. For small resolutions K this is most 
conveniently done by first computing the eigenvalues and then the matrix exponential. At larger 



resolutions we employ a random vector routine fll4fl to compute the trace of the matrix exponential, 



which has been approximated by Trotter decomposition. 

As a test case and to fix the range of external parameters T, L where reliable numerical results 
can be extracted we investigate the free Fermi gas. The light cone expression of the thermodynam- 
ical potential (density) of the free quantum gases is given as (upper sign fermion (/), lower sign 
boson (b)) 







2k V I V 2 2p 
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L/27I [MeV '] L/2JI [MeV"'] 

(a) whole range in L (b) small values of L 



Figure 1: The thermodynamical potential density — 2nT\n 3? as a function of L for different temperatures. 
In figure (a) the whole interval in L is shown while in figure (b) L is limited to small values to present the 
finite size effects. Open symbols in both figures are the potential at maximal resolution K = 110. Closed 
symbols are given by an extrapolation. The slope of the linear part (values selected colored in red) is fitted 
to extract the invariant potential density. 



Equation (2.5) is derived analogous to the instant form case, replacing the spatial volume by the 
light-like extension. In the large 'volume' limit the densities are equal in both relativistic forms. 
Figure l(a} shows results for the free electron gas of mass m = 0.5 eV at resolution K = 1 10. At 



small system volumes clear finite size effects are visible (see figure 1(b)) and at large volumes 



there are derivations from the exact result (2.5) because of the finite resolution. Therefore one has 



to identify a scaling window where the linear behavior in (2.4) shows up. Finding such a window is 
easy at small temperatures, but for increasing temperatures the scaling window is pushed to regions 



of large volumes. For the largest temperature shown in figure 1(a) the relative error is below 1.5%, 
see [ |l2| ] for more details. 



3. QEDi + i at finite Temperature on the Light Front 



The light front Hamiltonian of the massive, chiral Schwinger model (QEDi + i) is given in [ |11[ ] 
without the dynamical gauge field zero mode. Generically the Hamilton operator has the structure 



H 



m z H +—V 
71 



m z 1 



(3.1) 



where Hq is the free Hamiltonian, that is diagonal in free particle basis and V some complicated 
operator containing combinations of four creation and destruction operator of fermions and anti- 
fermions. The application of DLCQ to thermodynamics requires rather larger harmonic resolutions, 
as a byproduct one gets more accurate estimates for mass spectrum for different couplings. Still 
one has the extrapolate the raw data to the limit K —> °°, which was done by second-order power 
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K 1 / K 

(a) Full mass spectrum at m/g = 1 (b) Lowest bound states at m/g = 1 

Figure 2: The invariant mass spectrum at m/g = 1. Part (a) shows the full mass spectrum up to K = 35 and 
in (b) the six lowest mass eigenvalues of QEDj + i are depicted. The dashed line is a quadratic fit to the data 
(values used in the fit are colored in red) and used to extract the continuum limit. Mo is defined such that the 
continuum value of the lowest mass is normalized to one. 



functions in l/K. In figure || the mass spectrum is plotted for m/g = 1, which is in the non- 
perturbative coupling regime. Thereby 2(a) shows the full spectrum up to K = 35 and the growth 
of DLCQ states is apparent. For M/Mq > 2 the spectrum is continuous and we singled out the six 
lowest mass states in figure 2(b)| . 

A comparison with masses obtained by other means like finite lattice calculations [15], vari- 
ational DLCQ [ |l6| ] and fast moving frame approach [17] is possible for the lowest two states and 
our results [12] are generally in very good agreement. Slight differences appear for m/g < 2~ 3 . 
The reason is that the choice of the fermionic Fock representation is presently not optimal in this 
case. 

The thermodynamic quantities are obtained in the way outlined in the former section. Follow- 
ing two relations hold 



-ft) 



T 

-In St 



and // 



dlnSTT 2 
dT 77 



(3.2) 



In practice, we have directly computed the pressure and used a numerical derivative to determine 
the internal energy density u. Figure 3(a)| ( 3(b) ) show the dimensionless ratio of pressure (internal 
energy density) and T 2 as a function of temperature of a QED gas for four different couplings. 
For massive fermions we meet the chargeless condition of physical states Q\phy) = and thus 
computed 3f in the canonical ensemble. In the limit of vanishing mass we used grand canonical 
ensemble with jj. = since in this case LF QEDi + i is a free boson theory of mass mt, = -^=. The 
errors in figure 3(a) are due to the extrapolation to larger resolutions and roughly carried over from 
the free case computation in section 2. In figure 3(b' the fluctuations of the data points can be 
reduced by setting a smaller temperature grid. The external parameters T, L are both given in units 
of the lowest bound state mass Mq. Unlike to the free case before we do not set a definite physical 
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0.6 



0.2 




■ m/g = 2 

— ml° = Q 

■ > g/m = 



T/M 

(b) energy density 



Figure 3: The thermodynamical quantities pressure (a) and internal energy density (b) divided by T as 
functions of temperature. Four different couplings are shown: pure Bose gas (solid line), strongly interacting 
Fermi system m/g — 2~ 3 (squares), weaker interacting Fermi system m/g — 1 (circles), free Fermi system 
g/m = (diamonds). 



scale since Mq is not fixed in physical units. Remind that the mass of the first bound state can be 
large, like in QCD where the lightest bound state is the pion of mass m n = 140 MeV made out of 
nearly massless quarks. To judge whether the temperature reached in the numerical computation is 
sufficient we compare the high-temperature values in the figures || with the T S> regime of (|3|). 
One finds that the pressure has not yet reached the value expected by the high-temperature limit 
p/T 2 w 7T/6, but the internal energy is at T > Mq in the range of u/T 2 m n/6. In comparison to 
the results of the earlier study [^] p and u are computed at significantly higher temperatures and 
no sign of the conjectured phase transition is found. The figures || offer the interpretation that the 
thermodynamical quantities change smoothly under variation of the coupling. 



4. Conclusion 



This contribution is concerned with the application of light cone quantization to the thermo- 
dynamics of non-perturbative quantum field theory. As an example we treated QED in 1 + 1 di- 
mensions and presented the pressure and the internal energy. Since we have computed the partition 
function other thermodynamical quantities like entropy and the specific heat can be obtained via 
derivatives of In i2° and the equation of state can be given numerically. This procedure is limited 
by the exponential growth of basis states and the dimensionality of the Hamiltonian matrix with in- 
creasing harmonic resolution. To this end an effective (non-perturbative) renormalization program 
for Hamiltonians is necessary. Promising suggestions to this direction are the similarity transfor- 



mation renormalization transformation Jl9| , [20p, and the density matrix renormalization group in 
momentum space [^l]]. More specifically within the massive Schwinger model an inclusion of the 
dynamical zero mode is desirable because the condensate connected to the zero mode may have 



6 



Thermodynamical properties ofQED in 1+1 dimensions within light front dynamics 



Stefan STRAUSS 



impact on the thermodynamics. In Ref. [22] such a light cone Hamiltonian is suggested which 
could be a good starting point. 

A main objective of this direction of research is the extension to four-dimensional finite density 
QCD, avoiding the Monte-Carlo sign problem and reveal how the phase diagram of quark matter 
is from first principles. So far we have proven that the consistent application of the theoretical 



framework outlined in [ |13| ] to the non-perturbative situation is possible and leads to reasonable 
thermodynamic results. The numerical issues faced will surely increase, if one considers the full 
3+1 case. Nevertheless, the chance of eventually arriving at results for the phase diagram of QCD 
alternative to the ones of the well established lattice QCD is really exciting and along the demands 



recently claimed by Ken Wilson [23]. 
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